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We investigate tensor mesons as quark-antiquark bound states in a fully covariant Bethe-Salpeter 
equation. As a first concrete step we report results for masses of J PC = 2 ++ mesons from the chiral 
limit up to bottomonium and sketch a comparison to experimental data. All covariant structures of 
the fermion-antifermion system are taken into account and their roles and importance discussed in 
two different bases. We also present the general construction principle for covariant Bethe-Salpeter 
amplitudes of mesons with any spin and find eight covariant structures for any J > 0. 
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I. INTRODUCTION 



In QCD, mesons are viewed as bound states of 
(anti)quarks and gluons. Starting with a qq-picture they 
appear simpler than baryons and thus represent prime 
targets for theoretical investigations. Spin and the cor- 
responding meson degrees of freedom are essential for an 
understanding of the meson spectrum and properties in 
general. 

In a constituent-quark model (e. g. [IHl]), mesons with 
total spin J are easily obtained via adding units of or- 
bital angular momentum to a quark-antiquark state. In 
particular, given the quantum numbers J PC for a meson 
with equal-mass constituents, the parity P is given by 
(— and the C-parity C by (— Furthermore, 
the total spin J, the internal (quark-antiquark) spin s, 
and the orbital angular momentum I and their projec- 
tions have to satisfy the well-known addition rules for 
angular momenta. 



In the context of the Bethe-Salpeter equation (BSE), 
the Lorentz covariant structure of meson amplitudes 
(also for arbitrary spin) has in the past been investi- 
gated mainly in setups involving reductions of the BSE 
(e.g. )■ Herein we present the first covariant study 

of tensor mesons that is consistent with respect to the 
axial- vector Ward-Takahashi identity in the context of a 
Dyson-Schwinger-Bethe-Salpeter approach to QCD. 

The paper is organized as follows: Sec. [TTI sketches the 
formalism used and the corresponding details of immedi- 
ate necessity, Sec. lIIII contains the explicit construction of 
the covariant amplitude for a 2 ++ meson, the construc- 
tion principle for J > 2 amplitudes is given in Sec. IIV1 
the 2 ++ results are presented and discussed in Sec. W\ 
and we conclude in Sec. EU All calculations have been 
performed in Euclidean momentum space. 
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II. MESONS FROM THE BSE 

In this work, we employ (DCD's Dyson-Schwinger- 
equations (DSEs) (see, e.g. [2, [01 for recent reviews) 
together with the quark-antiquark Bethe-Salpeter equa- 
tion (BSE). The latter is the covariant bound-state equa- 
tion for the study of mesons in this context [7j . An analo- 
gous covariant approach to baryons is possible in a quark- 
diquark picture (e. g. I2 CH22J and references therein) or a 
three-quark setup |23l. |24||. 

While the goal of a self-consistent solution of all DSEs 
can be held up in investigations of certain aspects of the 
theory (see, e.g. [2ol |26| and references therein), numer- 
ical hadron studies such as ours require employment of a 
truncation. For our first covariant look at tensor mesons 
we use the so-called rainbow-ladder (RL) truncation. It 
is both simple and offers the possibility for sophisticated 
model studies of QCD within the DSE-BSE context, 
since it satisfies the relevant Ward-Takahashi identities 
(WTIs), namely the axial- vector (see e.g. [23, l28j ) and 
vector (see e.g. [2!|[3(|) WTIs. The literature regarding 
the employment of terms beyond RL truncation can be 
traced back from e.g. [U, 32] . The axial- vector WTI is es- 
sential to see chiral symmetry and its dynamical breaking 
correctly realized in the model calculation from the very 
beginning. As the most prominent result, one satisfies 
Goldstone's theorem [23] and obtains a generalized Gell- 
Mann-Oakes-Renner relation valid for all pseudoscalar 
mesons and all current-quark masses [33l |3J) . We note 
that this relation can be checked numerically and is sat- 
isfied at the per-mill level in our calculations. 

The general structure of the BSE for a meson with spin 
J, total qq momentum P and relative qq momentum k 
or q, respectively, is 

r^ -(k;P) = [ K(k;q;P)S(q + )T^(q;P)S(q^ ,(1) 

Jq 

where the semicolon separates four- vector arguments. 
Y» u -(k;P) is the Bethe-Salpeter amplitude (BSA) and 
has J open Lorentz indices \w . . .. The dressed-quark 
propagator S(p) is obtained from the quark DSE, the 
QCD gap equation. Since our focus here is the BSA, 
we refer the reader to [33, uia KSgl for more details on 
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the quark DSE and to [13] for a description of our 
corresponding numerical solution method. In the BSE 
the quark and antiquark propagators depend on the 
(anti)quark momenta q + = q + nP and g_ = q— (1 —rj)P, 
where n £ [0, 1] is a momentum partitioning parameter 
usually set to 1/2 for systems of equal-mass constituents 
(which we do as well), J — j A d i q/(2n) i represents 
a translationally invariant regularization of the integral, 
with the regularization scale A [33[ . 

The kernel K in the homogeneous, ladder-truncated qq 
BSE is essentially characterized by an effective interac- 
tion G(s), s := (k — q) 2 . Following [36|, an ansatz used 
extensively for many years [35| is employed here, which 
reads 

G(S) _ 4lT 2 D /u) 2 47T lm 7TT( S ) 

l/21n[r+(l + S /A| CD )2]- ^ ' 

This form provides the correct amount of dynamical 
chiral symmetry breaking as well as quark confinement 
via the absence of a Lehmann representation for the 
dressed quark propagator. Furthermore, it produces 
the correct perturbative limit, i. e. it preserves the one- 
loop renormalization group behavior of QCD for so- 
lutions of the quark DSE. As given in [35], J-(s) — 
[l-exp(-s/[4m 2 ])]/s, m t = 0.5 GeV, r = e 2 -l, N f = 4, 

Aq£= 4 = 0.234 GeV, and 7m = 12/(33 -2JV». Note that 
the same effective interaction appears also in the corre- 
sponding rainbow-truncated quark DSE. 

This function, which mimics the behavior of the prod- 
uct of quark-gluon vertex and gluon propagator, is 
mainly phenomenologically motivated. While currently 
debated on principle grounds (e.g. |26|,[3j|) the impact of 
its particular form in the far IR on meson masses is ex- 
pected to be small (see also (39] for an exploratory study 
in this direction). 

D and w, in principle free parameters of the model in- 
teraction, can be used to investigate certain aspects of 
both the interaction and the bound states in the BSE. In 
particular one can interpret D as an overall strength and 
lu as an inverse effective range of the interaction (for more 
details and a thorough discussion of parameter depen- 
dence of the results see (3(|), a notion first investigated 
in the study of radial meson excitations [U, H(| ■ In the 
range ui £ [0.3, 0.5] GeV, the prescription D xuj = const. 
follows from fitting of the model parameters to ground- 
state properties [35J and defines a one- parameter model, 
which is the setup and range used in [36| and also here. 
With all ingredients specified, the BSE is solved numer- 
ically, a procedure well under control |4lj . 

III. TENSOR-MESON BETHE-SALPETER 
AMPLITUDE 

The BSA r M "'"(fc; P) of a meson as a bound state of 
a quark- antiquark pair depends on two four-vector vari- 
ables: the total as well as the relative quark-antiquark 



four-momenta P and q, respectively. They can be param- 
eterized in terms of the Lorentz-invariant scalar products 
P 2 , q 2 , and q ■ P. The fermion-antifcrmion spin proper- 
ties are encoded in the 4x4 matrix structure of Y^- - 
[7j , where the open Lorentz indices appear in connection 
with the total spin of the state. A corresponding basis 
of linearly independent structures {T^ v "} (i — 1, . . . , N) 
involving Dirac matrices allows one to expand the BSA 
into a sum of Dirac covariants and the corresponding 
scalar coefficients Fi , which we will subsequently refer to 
as components [4l| . The latter only depend on the afore- 
mentioned scalar products P 2 , q 2 , and q ■ P, and one 
gets 

N 

r^-(fc;P; 7 ) =£2f-(fc;P; 7 ) F(q 2 ,q-P,P 2 ) ,(3) 
»=1 

where the dependence on 7" has been made explicit and 
a generalized scalar product for the covariants Tjf "' is 
defined via the Dirac trace 

X lr / - i " 7 T" '-,/:'-./*• (4) 

. . . 

One may also choose the basis elements orthogonal such 
that tij = Sij, with the f(i,j) functions of q 2 , P 2 , and 
q ■ P, or orthonormal such that in addition f(i, j) = 1 for 
all i, j. The sum is carried out over the J indices u, v, . . .. 

Note that for an on-shell BSA P 2 = -M 2 is fixed, 
while one artificially varies P 2 in the solution process of 
the homogeneous BSE. In the corresponding inhomoge- 
neous BSE one has P and therefore also P 2 as a com- 
pletely independent variable (see, e.g. (4lTl43l |). 

Thus, the on-shell scalar components Fi(q 2 ,q ■ P,P 2 ) 
effectively depend on the two variables q 2 and q ■ P, the 
latter of which can be parameterized by the variable z G 
[—1,1] related to the cosine defining the angle between 
the four vectors P and q. In principle, the components 
Fi can be expanded further in Chebyshev polynomials, 
but we do not use such an expansion here (for details and 
an illustration of Chebyshev moments, see [H, H3|). 

With the independent four momenta and 7" one can 
construct four independent Lorentz-scalar structures, 

1, 7-P, 7-9, P , (5) 

where cr q,p :— i/2 [7 • q, 7 • P]. These four covariants, 
which provide a basis corresponding to scalar mesons 
( J p = + ), serve as the basic building blocks for any me- 
son BSA. Together with pseudoscalar covariants ( J p = 
0~) as well as the bases for J = 1 for all corresponding 
quantum numbers, these were explicitly constructed in 
[3a |. Here we concentrate on J = 2 and higher. 

For J = 2 one has 8 independent covariant structures 
in the BSA. Let 

% ■= % -P^-jyT ' ( 6 ) 

ll -=^-P^- > (7) 

tt 7 ' P r 7 ' Q T 

7^ : =7m-^u^2--^77tW > ( 8 ) 
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be transverse projections of 7 and q with respect to the 
total meson momentum P and each other (in particu- 
lar the vectors {P M , q J , 7j T } are orthogonal to each 
other). 

Defining furthermore the transverse projection of the 
metric 



T r _ PyPv 

9pu — °nv p 2 



(9) 



and the two transverse, symmetric, and traceless struc- 
tures 



= + 9J7J - ~gl„ 7 • q T (10) 



(11) 



one obtains the following set of tensor ( J p = 2 + ) covari- 
ants 0] 



T^" = AP" j-qq-P- 2N^ q ■ P 



r['" = .M 1 "- , , 
W" 7 • P 



rpflV 

J 3 — 



T% v = 2M» V a q P - UN"" 7 • P 



1 7 — 



i 7 • P q ■ P 



1 a — 



-2iN^ 



(12) 



Note that T 5 . . . T 8 were only given implicitly in Q • All 
Ti as given here are even under charge conjugation (for 
details, see e.g., (Hill). Thus, to obtain a J PC = 2++ 
state, all components Ft must be even functions of q ■ P, 
which for the present setup is indeed the property of the 
ground state in the system. Note also that these covari- 
ants are in general neither orthogonal nor normalized; 
orthonormal covariants can be generated via a Gram- 
Schmidt procedure applied to the set of terms in ([5]>. 
leading to 



7 • P, j-q 



ri.P 



(13) 



To orthogonalize the above 2 + covariants one introduces 
the symmetric and transverse expressions 



M, 



in 

= q^ql 



7j T <?J + <?j7j T and 



(14) 
(15) 



which automatically satisfy Eq. The next step is 

to implement the tracelessness, which is equivalent to 
orthogonality with respect to g£ v . This yields 



= M M „ - g 



-p Mpa9pa 
p N pa 9 pa 



(16) 
(17) 



which corresponds to Eqs. (fTO)) and ([TT|) . and by mul- 
tiplication with the four scalar covariants in (fTB"]) gives 
the eight desired orthogonal tensor covariants. Note, 
however, that Eqs. (fTU|) and ([To]) are slightly differ- 
ent. Subsequently, normalization is achieved via X!j = 



Ti/y/Tr[Ti-T^ 




FIG. 1. (Color online) dependence of meson masses on m w 
(the pseudoscalar-meson mass calculated for a given current- 
quark mass). Vertical dotted lines correspond to positions for 
light, strange, and charm qq states. 



IV. BSA FOR ANY MESON SPIN 

To consider mesons of any particular spin J, one has 
to construct Lorentz-tensors of rank J which are totally 
symmetric, transverse in all open indices and Lorentz- 
traceless (see, e.g., 45]): such an object has the 2 J + 1 
spin degrees of freedom as demanded in quantum me- 
chanics of a massive particle. These restrictions, together 
with the properties of the Dirac matrices, lead to eight 
covariant structures for J > 1. More precisely, the two 
tensors M^ v and defined above can be generalized 
such that Nfj, v ,,, T is the traceless part of 



T T 
% 1v 



T 
<1t 



(18) 



and M Mt/ ... T is the traceless part of the totally symmetric 
sum constructed from 



TT T T 
Ip 1v ■■■ Qr 



(19) 



Each of these multiplied by the four terms in ([13]) defines 
four rank-J tensor covariants, in total eight, orthogonal 
in the sense of Eq. 

Obviously, Eqs. (fT5|) and (fT7|) follow from this con- 
struction. As a further quick check we consider the sim- 
plest such example, namely a vector meson: from J = 1 
one immediately obtains = <?J to give the first four, 
and M M = 7J 71 to give the second four covariants. 



V. RESULTS 

Here we present results for J PC = 2 ++ states that 
extend the study of Ref . [36[ . Consequently, we present 
correspondingly augmented figures here. Fig.[T]shows the 
meson masses for pseudoscalar, scalar, vector, axialvec- 
tor, and tensor qq states as functions of the pion mass, 
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u/d s c b 




0.3 0.4 0.5 0.3 0.4 0.5 0.3 0.4 0.5 0.3 0.4 0.5 

CO CO 0) CO 

FIG. 2. (Color online) dependence of meson masses on cj. 
Dotted lines correspond to experimental data [46l ] where ap- 
plicable. Figure updated from Ref. [36j 



obtained from the BSE in RL truncation employing the 
effective interaction of Eq. @. The three vertical dotted 
lines indicate the positions of the nn, ss, and cc states, 
respectively (in the usual notation, n here denotes light 
quarks). Note that in RL truncation one cannot eas- 
ily employ arbitrary flavor mixing between 5(7(3)-flavor 
octet and singlet states; all states are either purely nn or 
ss, which corresponds to ideal mixing. As expected, the 
2 ++ mass lies above all other states for the entire range 
from the chiral limit to bottomonium. 

Comparison to experimental data is shown in Fig. [5] 
via the dotted lines in appropriate colors, where for the 
nfi, ss, cc, and bb cases separately the dependence of 
the bound-state masses on the model parameters is stud- 
ied. More precisely, as mentioned in Sec. [Ill the inverse 
effective range ui is used to explore the state's sensitiv- 
ity to the details of the long-range part of the strong 
interaction [34], HH- Two observations are noteworthy: 
First, the 2 ++ mass shows the same w-dependence as the 
other orbital excitations for each of the four columns. 
Secondly, the agreement with experimental data is sig- 
nificantly better than for the often-quoted axial-vector 
states. In the case of the latter, reconciliation of an RL 
study constrained by pseudoscalar- and vector-meson ob- 
servables seems unlikely, while this need apparently not 
be the case for tensor mesons, indicating that the lat- 
ter are simpler and not as sensitive to the details of the 
quark-gluon vertex employed in a DSE-BSE study as the 
axial- vector states. 

A further technical note concerns the 2 ++ results for 
u) = 0.5 GeV: Due to the analytic structure of the quark 
propagators for this parameter choice, the masses of the 
2++ mesons are only accessible to us via extrapolation 
techniques (see [47] for a discussion). However, the ex- 



trapolations used are reliable and stable, and the result- 
ing error bars are smaller than the size of the symbols in 
Fig. [5] except for the u/d case, where we get an uncer- 
tainty of ±75 MeV. 

An interesting question related to the issue of "sim- 
plicity" of the 2 ++ states in this approach is, how impor- 
tant the various covariants/components are in the BSA 
or, in other words, how many covariants are needed to 
arrive close to the full result. One possibility to inves- 
tigate this is to leave out each covariant and recompute 
the mass of the state with the remaining seven. Small 
differences to the full result then indicate covariants of 
minor importance. Naturally there is a caveat for such 
an investigation, namely that the choice of the covariants 
is somewhat arbitrary. 

In our case we used two sets of covariants: the one 
given explicitly above in Eq. (|12[) and the other, orhonor- 
mal, constructed according to the principles detailed in 
Sec. IIVI We have performed this test for both sets of 
covariants and present the results in Tab. fl] We enumer- 
ate the orthonormal covariants in the following way: the 
four terms in (fT3")l multiplied with (IT71) are numbered 1 
to 4, and (fT3"|) multiplied with (| 16[) yield covariants 5 to 
8. For either set, one needs five of the eight covariants 
to arrive at a number which is within one percent of the 
full result. Furthermore, omitting the contribution from 
N^ v as indicated in |7| for this particular case yields a 
number which is 7% too low compared to the full result. 

VI. CONCLUSIONS AND OUTLOOK 

We have presented the complete set of Dirac covariants 
for mesons of spin 2 and given an explicit construction 
principle for the corresponding set of covariants for an 
arbitrary spin J for the first time. We have furthermore 
explored 2 ++ states in a well-established RL truncated 
model setup of QCD's DSEs and solved the correspond- 
ing quark- antiquark BSE numerically. The results are 
both reasonable and surprising in that they follow ex- 
pected patterns, but are closer to experimental data than 
axial vector mesons, even in the present simple setup. 

The numerical calculation of further states with J PC = 
2 3 , etc. is work in progress and will be presented 
in future publications. Naturally, this includes radial ex- 
citations of these states and opens up the concrete pos- 
sibility to investigate Regge trajectories in the covariant 
DSE-BSE approach. 
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TABLE I. Meson mass of the ss 2 ++ state with ui = 0.4 GeV with all covariants included as well as with single covariants left 
out. The change in bound-state mass is given compared to the full result. The results are presented for both the covariants of 
Eq. (|12[) and the orthonormal set of covariants constructed thereafter. All numbers are given in GeV. 



Covariant 


missing 


none 


1 


2 


3 


4 


5 


6 


7 


8 


Eq. CE2) 


Mass 


1.448 


1.575 


1.455 


1.502 


1.509 


1.502 


1.287 


1.452 


1.450 




Change 


+0.000 


+0.127 


+0.007 


+0.054 


+0.061 


+0.054 


-0.161 


+0.004 


+0.002 


Orthonormal 


Mass 


1.448 


1.502 


1.445 


1.540 


1.420 


1.669 


1.457 


1.446 


1.508 




Change 


+0.000 


+0.054 


-0.003 


+0.092 


-0.028 


+0.221 


+0.009 


-0.002 


+0.060 
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